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Abstract 

The time evolution of the entanglement of a pair of two spin qubits is investigated when the 
two qubits simultaneously couple to an environment of an anisotropic Heisenberg XX Z spin chain. 
The entanglement of the two spin qubits can be created and is a periodic function of the time if 



^ the magnetic field is greater than a critical value. If the two spin qubits are in the Bell state, the 

Oh 

+!> entanglement can be stored with relatively large value even when the magnetic field is large. 
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I. INTRODUCTION 



Entangled quantum states are used mainly for quantum information processing, such as 
quantum teleportation, quantum secret-code and quantum computation [IH3] . Many investi- 
gations showed that entanglement exists naturally in the spin system when the temperature 
of the system is at zero [IH7]. In recent years, the study of the dynamics of entanglement 
[EHI5] has attracted much attention in the manipulation of quantum systems. The dynami- 
cal properties and the time evolution of the entanglement in different quantum systems were 
investigated. These systems included mobile particle elastically-scattered by static spins [8], 
quantum mixed states [91 [10], two oscillators coupled to the same environment [11, 12J, two 
d-level systems [13J, decoherence of a spin-1/2 particle coupled to a spin bath in thermal 
equilibrium [14J, a spin chain in driving the the decoherence of a coupled quantum system 
[To] , etc. Meanwhile, the effects of the environment were taken into account. The excitation 
and quantum information transfer was investigated between two external spins when they 
coupled to a one-dimensional spin chain at different sites [16J. The entanglement induced by 
two external spins could be used to signal the critical points when they were simultaneously 
coupled to an environmental XY spin chain [TTl ITS] . The decay of the Loschmidt echo was 
enhanced by the quantum criticality of the surrounding Ising chain when an external spin 
was coupled to the environment [19] . When two external spins coupled to a transverse field 
Ising chain, the induced entanglement could be enhanced near quantum criticality and could 
be used to detect the quantum phase transition (20] , which occurred in the many-body quan- 
tum systems [21]. The dynamical properties of the entanglement in a spin system need to 
be further investigated when it is coupled to an external environment. 

In this paper, the dynamics and the time evolution of the entanglement of a pair of two 
qubits are investigated when the two qubits simultaneously couple to an environment of an 
anisotropic antiferromagnetic Heisenberg spin chain with magnetic field. In Section II, the 
Hamiltonian of the system and the effective Hamiltonian of the two qubits coupled to the 
environment are presented. In Section III, the time evolution of the system is analyzed for 
the simplest case of the environment. The entanglement creation in the coupled pair of two 
external qubits is discussed in Section IV. In Section V, the storage of the entanglement in 
the coupled pair of two external qubits is investigated. A discussion concludes the paper. 
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II. HAMILTONIAN OF THE SYSTEM 



When two external spin qubits are coupled with the environment of a one-dimensional 
spin chain, the Hamiltonian of this system can be written as 

H = H + Hj. (1) 

where H is the Hamiltonian of the environment. If the environment is an anisotropic Heisen- 
berg XX Z spin chain, one has 

N N 

Ho = J + °Mu + <f i) + B E > ( 2 ) 

i=l i=l 

were cr"(a = x, z) are Pauli operators, N is the number of the spin chain, J is the coupling 
coefficient between the spins, B is the magnetic field along the z-axis with the anisotropy 
Aj = A G (0,1). In Eq. (1), Hi is the interaction Hamiltonian between the two external spin 
qubits and the environment and can be written as, 

N 

Hi = J p ^(OaCTi + a b ai). (3) 
1=1 

where a a and o\, are the Pauli operators of the qubits a and b, J p is the coupling coefficient 
between the external spin qubits (a and b) and the Heisenberg spin chain. In order to 
facilitate the calculation, the coupling coefficients are chosen as J = 1 and J p = 0.2 in this 
paper. That is, the environment is represented by the antiferromagnetic Heisenberg XXZ 
spin model. The schematic diagram of the system is shown in Fig. 1. The two qubits are 
symmetrically located at the two sides of the spin chain. 

Frohlich transformation [201 122] can be used to solve the problem of induced effective 
interaction between two qubits through the medium of the Heisenberg spin chain. The 
environment of the antiferromagnetic Heisenberg spin chain has non-degenerate ground state 
|^o) with ground state energy i?o- According to the standard canonical transformation [201 
1221 123] , the effective Hamiltonian of the external spin qubits can be written as 

E 1l = Xl E S -E • (4) 
where the projector is Pj = \i>j){^j\ and = 1,2,... A;) is the time dependent excited 

state with energy Ej. After some straightforward calculations, the effective Hamiltonian can 
be reduced to 
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h: 



eff ~ 



ab 




where the parameters are m, 




the real part of the product (m Q n*) with a,(3 = x,y,z. When the eigenstate \ipj) and the 
corresponding eigenvalue Ej of H are obtained, the effective Hamiltonian H^tf can be easily 
calculated. 

III. ANALYSIS OF TIME EVOLUTION 

In order to describe the time evolution of the entanglement of two-qubit system, the 
concurrence is used as a measure of the entanglement. The concurrence is defined as [24] ES] 



where the Aj(i = 1, 2, 3, 4) are the square roots of the eigenvalues of the density matrix g ab . 
The density matrix g ab is given by 



The ground state of the environment of the Heisenberg spin chain can be chosen as |</>o) 
while that of the two external spin qubits a and b can be chosen as |01). Under the influence 
of the environment, the two external spin qubits have an initial state as follows 



C = max{Ai — A 2 — A 3 — A 4 , 0}, 



(6) 



Qab = Pn(crl ® o-f)pi 2 (af ® o%). 



(7) 



|^ ) = |0 O ) (g, |01) a6 . 
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The time evolution of the state is 



m))=exp(-zH: b ff t)\iP Q ) ab , 



(9) 



with the density matrix g ab = \if)(t))(if)(t)\. 
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The reduced density matrix g a b(t) can be written as 

/ u(t) \ 

Qab{t) = (10) 

, y*(t) w 2 (t) 
\ v(t) J 

in the standard basis {|00), |01), |10), |11)}. The corresponding concurrence C(t) of the two 
external spin qubits can be calculated from the reduced density matrix g a b(t) and given by 



C(t) = 2max{|y(t)| - y/u(t)v(t),0}. 



(11) 



IV. ENTANGLEMENT CREATION 



For the simplest case of N = 2 in the anisotropic Heisenberg XXZ spin chain, the 
eigenenergies and eigenstates of the system are E\ — A — 2B,E 2 = A + 2B,E 3 = —A + 
2,E 4 = -A -2 and = = |00)> 3 ) = ^(|01) + |10)), \<p 4 ) = ^(-|01) + 

1 10) ) respectively When B — A > 1, the ground state is |0 O ) = Then the effective 

Hamiltonian can ^ e wr itten as 



Tjab 



where the parameter g is given by gf 



( I o\ 

12 

2 10 

\0 1 / 

J; 



P (B-A+1)(B-A-1) ' 



(12) 



The matrix elements 



u(t),w 1 (t),y(t),w 2 {t) and v(t) in Eq. (10) are given by u(t) = 0,wi(t) = \ + le"* 4 ^ + 



l e ^,y(i) = -\e-^ + le i4 ^, W2 (t) = 



l p -i4gt 
4 C 



|e i49 *,w(t) = 0. When B - A < 1, 



the effective Hamiltonian -f^?// = 0- There is no entanglement between spin qubits a and b. 
When B — A = 1, the ground state energy equals the excited state energy, i. e., E\ = E 4 . 
The energies of the two states are crossed at this point. Since the two states are degenerate, 
Eq. (4) is not valid to calculate the effective Hamiltonian -ff"// wnen B — A — 1. That is, 
there is a critical value of the magnetic field B c . The value of B c is given by B c = 1 + A. 
If B < B C) the concurrence C(i) is zero. That is, there is no entanglement when B < B c . If 
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B > Bq, the entanglement appears. That is, the entanglement can be created when B > Be- 
Then the concurrence C(t) can be given by 



The concurrence C(t) as a function of the time t is plotted in Fig. 2 when the magnetic 
field B and the anisotropy A are varied. The values of the anisotropy are A = 0.2, 0.4 and 
0.6 with B > Be in Figs. 2(a), 2(b) and 2(c) respectively. From Fig. 2, it is seen that the 
concurrence C(t) is a periodic function of time t. It almost oscillates between the maximum 
value of one and the minimum value of zero. The period decreases as the magnetic field B 
increases. 

The anisotropic antiferromagnetic Heisenberg XXZ model was used to investigate the 
order-to-disorder transition of the material Cs^CoCl^ [26]. For the material CS2C0CI4, the 
anisotropy is A = 0.25. When the number of spins in the environment of the Heisenberg 
XXZ chain is greater than two, there is no approximate analytic solution of Hflf and 
C(t). To calculate C(t), the numerical computation needs to be performed. In Fig. 3, the 
concurrence C (t) is plotted as a function of time t when the spin numbers in the environment 
are N = 4, 6 and 8. From Fig. 3, it is seen that the concurrence C(t) is a periodic function 
of t with two different kinds of periods. Both periods decrease as the spin number N in the 
environment increase. There is a critical value Be of the magnetic field. When B < Be, 
the concurrence C(t) oscillates following the large period. The period decreases slightly as 
B increases. While B > Be, C(i) oscillates following the small period. The period increases 
as B increases. The concurrence C(t) can be approximately given by 



Where g(N) is a function of the spin number N in the environment. Though there is no 
analytic expression of the critical field Be, it can be numerical calculated. The critical field 
Be is plotted in Fig. 3(d) as a function of 1/N. From Fig. 3(d), it is seen that Be decreases 
linearly as 1/N decreases. In the thermodynamic limit of N — > 00, Be tends to zero. The 
regime for larger period of oscillation disappears. 




(13) 




(14) 



6 



V. ENTANGLEMENT STORAGE 



The concurrence C(t) is plotted as a function of magnetic field B and time t in Fig. 4 when 
the initial state of the two external spin qubits a and b is in the Bell state ^(|01) + |10)). 
The anisotropy is chosen as A = 0.25 [26]. The spin numbers in the environment of the 
anisotropic antiferromagnetic Heisenberg chain are N = 2,4, 6, and 8. From Fig. 4, it is seen 
that the concurrence C(t) is a oscillation function of time t. The oscillation period decreases 
as B increases. Obviously, the concurrence C(t) is divided into several regions by different 
critical values of the magnetic field Be- The red circles in Fig. 4 show the critical values of 
Be- At the critical point Be, the energies of the ground and excited states are crossing and 
the states are degenerate. For N = 2, there are two parts in C(t) divided by one Be- For 
B < Be, C(t) is almost a constant of C(t) = 1.0. For B > Be, C(i) oscillates with a small 
period [cf. Fig. 4(a)]. For N = 4, there are three parts divided by two different values of Be 
[(marked by two red red circles in Fig. 4(b)]. In the first part, C(t) is very close to 1.0. It 
oscillates with quite small amplitude. In the second part, C(t) oscillates with small period. 
In the third part, C(t) oscillates with even smaller period [cf. Fig. 4(b)]. When N = 6 
and 8 in Figs. 4(c) and 4(d), similar phenomena occurs. Obviously, the concurrence C(t) 
is divided into {N/2 + 1) parts by N/2 critical values of Be- The energy is crossing at the 
critical values of Be in the ground state as well as in excited states. The concurrence C(t) is 
jumping as the state is switched from an entangled state to another. In the thermodynamic 
limit, the continuous energy level crossings occur [27] . The first part of concurrence C(t) 
disappears. Other parts of C(t) tends to smooth and continuous. In Fig. 3(d), only the first 
critical value of Be as a function of spin number 1/N is plotted. From Fig. 4, it is also clear 
that the entanglement C(t) can keep large value even for relatively large magnetic field B. 
If the initial state of the two external spin qubits is the Bell state -^=(|00) + |11)), similar 
results as that shown in Fig. 4 are obtained. 



VI. DISCUSSION 



The time evolution of the entanglement of two external spin qubits is investigated when 
they are coupled to the environment of an anisotropic antiferromagnetic Heisenberg XX Z 
spin chain with magnetic field. The approximate form of the effective Hamiltonian is derived. 
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The concurrence is used as a measure of the entanglement. When there are two spins in the 
environment, there is no entanglement between two external spin qubits when the magnetic 
field is smaller than a critical value. When the magnetic field is greater than the critical value, 
the entanglement can be created and is a periodic function of the time. The entanglement 
almost oscillates between one and zero. The oscillation period decrease as the anisotropy 
and the magnetic field increase. There are N/2 + 1 parts in the entanglement divided by 
N/2 values of critical magnetic fields. The first critical magnetic field tends to zero when the 
spin number in the environment tends to infinity. When the initial state of the two external 
spin qubits is in one of the Bell state, the entanglement can be stored. Though there are 
different regimes in the entanglement, the entanglement always keeps quite large value when 
it oscillates with increasing number of spins in the environment. 
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Figure Captions 



Fig. 1 

The schematic diagram of two external spin qubits symmetrically coupled to the environ- 
ment of an anisotropic Heisenberg spin chain. 
Fig. 2 

The concurrence C(t) is plotted as a functions of the time t for N = 2 when the magnetic 
field B and the anisotropy A = are varied with B > Be- (a). A = 0.2. (b). A = 0.4. (c). 
A = 0.8. 

Fig. 3 

The concurrence C(t) is plotted as a function of the magnetic field B and the time t for 
different spin numbers N of the environment in (a), (b), and (c). The anisotropy is A = 0.25. 
(a). N = 4. (b). N = 6. (c). N — 8. (d). The critical field B c is plotted as a function of 
1/N. 

Fig. 4 

The concurrence C is plotted as a function of the magnetic field B and the time t for 
different spin numbers in the environment. The anisotropy is A = 0.25 and the Bell state 
-W|01 > +|10}) is chosen, (a). N = 2. (b). N — A. (c). N = 6. (d). N = 8. 
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Fig. 1 
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B=1.4 

(a)A=0.2 b=i.6 (b) A-0.4 
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Fig. 2 
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Fig. 4 
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